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ABSTRACT 

PSR J1734— 3333 is a radio pulsar rotating with a period P = 1.17 s and 
slowing down with a period derivative P = 2.28 x 10~ 12 , the third largest among 
rotation-powered pulsars. These properties are midway between those of normal 
rotation-powered pulsars and magnetars, two populations of neutron stars that 
are notably different in their emission properties. Here we report on the mea- 
surement of the second period derivative of the rotation of PSR J1734— 3333 and 
calculate a braking index n = 0.9 ± 0.2. This value is well below 3, the value 
expected for an electromagnetic braking due to a constant magnetic dipole, and 
indicates that this pulsar may soon have the rotational properties of a magnetar. 
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While there are several mechanisms which could lead to such a low braking index, 
we discuss this observation, together with the properties exhibited by some other 
high-P rotation-powered pulsars, and interpret it as evidence of a possible evo- 
lutionary route for magnetars through a radio-pulsar phase, supporting a unified 
description of the two classes of object. 

Subject headings: pulsars: general — pulsars: individual (PSR J1734— 3333) - 
stars: neutron 



Introduction 



Radio pulsars and magnetars are believed to be neutron sta rs which have been formed in 

the c ollapse of the cores of massive stars in supernova explosions (|Pacinilll967l ; lDuncan fe Thompson 
19921 ) . Despite this common origin, they have very different rotational and radiative prop- 
erties ( IWoods fe Thompson! l2006t iMereghettil 120081 ). While radio pulsars generally rotate 
with periods P of 0.01-1 second and exhibit no remarkable X-ray emission, magnetars have 
periods ranging from 2-12 seconds and are generally intense sources of X-rays, sporadically 
undergoing outburst episodes. 



The original theory of magnetars (IDuncan fe Thompson! 119921 ) suggested that some 
neutron stars evolve under different circumstances to normal, or rotation-powered, radio 
pulsars. Convection and fast rotation in their early stages would quickly generate strong 
magnetic fields, ultimately producing Soft Gamma Repe aters (SGRs) and Anomalo us X- 
ray pulsars (AXPs), the two manifestations of magnetars ( IWoods fe Thompson! 120061 ) . De- 
cay of their extremely-high magnetic field is believed to be the key to explaining outburst 
episodes and the excess of X-ray luminosity over the i r rate of loss of rotational energy 
( iThompson fe Duncan! Il995t IWoods fe Thompson! l2006t IMereghettil 120081 ) . 



However, it is not clear whether they are indeed born wit h different properties or whether 



2005; 


Gavriil et al. 


2008; 


Hales et al. 


2009) 



properties approaching those of magnetars and the existen ce of objects of both types having 



some similar emission properties is now w ell established (IKramer et al.l 120071 : iGavriil et al. 



2008; 


Camilo et al. 


2008; 


Levin et al. 


2010) 



of these intermediate objects, PSR J1734— 3333, in order to study its rotational evolution. 
In this letter we present the measurement of the braking index of this pulsar, a most unusual 
value which defies the standard models for pulsar spin-down. 
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1.1. Braking index and movement on the P — P diagram 



The rotational relationship between radio pulsars and magnetars is best demonstrated 
in the period - period-time-derivative (P — P) diagram (Fig. [1]). Radio pulsars are believed 
to be formed in the upper left region of the diagram, with short period and large slow-down 
rate, and move generally to the right and downwards. Magnetars have longer periods and 
larger slow-down rates than young radio pulsars and are found in the upper right region 
of the diagram. The chief rotational slow-down mechanism is often assumed to be loss 
of energy through electromagnetic radiation produced by a dipolar rotating neutron-star 
magnetic field. Thus, the position in the diagram is conventionally interpreted in terms 
of the constant surface dipolar magnetic field at the magnetic equator required to slow 
down a standard neutron star (B s = 3.2 x 10 G), as well as the chara cteristic age 

r c = —vj2i> = P/2P of the neutron star (with v = l/P)(e.g. lLyne fc Smithll2006l ). According 
to this, magnetars have surface magnetic fields of 10 14 ~ 15 G, about two orders of magnitude 
higher than those of most normal radio pulsars, and are as young as young radio pulsars 
(1-100 kyr). 

Motion in the P — P diagram can be described by the braking index, n, defined by 



v 



or 



P = kP 



2-n 



where k is a positive constant. With this definition, the slope of the evolutionary path in 
the diagram is equal to 2 — n. For pure magnetic braking with a constant dipolar magnetic 
field we expect n = 3. As a result, neutron stars are expected to follow lines of constant 
dipolar magnetic field (Fig. [I]), with a slope of —1. However, systematic variations of the mo- 
ment of inertia or magnetic field, non-dipolar braking resulting from higher-order magnetic 
multipoles or stellar winds could effectively change the slope of movement in the diagram 
jManchester et al.lll985l : iBlandford fc Romanilll988[ ). 



Observationally, the braking index can be determined from measurement of the first 
two derivatives of z/, since differentiation of the above power-law (Eq. [1]) gives n = vvjv 2 . 
Unfortunately, few values of n have been determined because young pulsars are often not 
stable rotators and short-term timing irregularities in the form of glitches and timing noise 
often mask any systematic motion in the P — P diagram. However, with long observational 
time-baselines it is often possible to establish the underlying movement across the diagram. 
So far, seven puls ars have steady e n ough rotation that stable values of braking index have 
been determined ( Lvne et al. 19931 . 1996 : Middleditch et al. 2006 : Livingstone et al. 2007 : 
Weltevrede et al.l 120111 ) (see Table [T]) . These values are all less than three and hence the 
slope of their movement in the P-P diagram is somewhat more positive than —1. The 
simple dipolar spin-down model is clearly imperfect. 



-4- 



2. The braking index of PSR J1734-3333 



PSR J1734— 3333 lies in the region between young pulsars and magnetars on the P — P 
diagram and is a young radi o pulsar (r r . = 8 kyr), p ossibly associated with the supernova 



remnant (SNR) G354.8— 0.8 ( [Manchester et al.ll2002l ). It has a high inferred dipolar mag- 
netic field (5 x 10 13 G), one of the highest known values amongst rotation-powered pulsars, 
and similar to the lowest valu es amongs t magn etars (even without considering the SGR 
with a very low P reported bv lRea et al.l (12010 ) ; see Fig. [I]). It was discovered during the 



Parkes Multibeam Pulsar Survey (IMorris et al. 



2002j), and it has been observed regularly 



since August 1997 by the 64-m telescope at Parkes and the 76-m Lovell telescope at Jodrell 
Bank. Like other high-P radio pulsars, PSR J1734— 333 3 presents modest X -ray luminosity, 
not surpassing its rotational spin-down energy loss rate (lOlausen et al.ll2010l ). although with 
somewhat higher temperature than is obs erved for pulsars with lower effective dipolar mag- 
netic fields but similar characteristic age (jZhu et al.ll2010l ). It exhibits timing noise typical 
of a pulsar of its age and has not glitched in 13.5 years, making possible a measurement of 
the v due to secular slow-down. 

The evolution of the timing residuals of PSR J1734— 3333 relative to a simple slow-down 
model of the rotational frequency and first derivative is shown in Fig. [2j The curve is a fit 
for a second derivative as well, and these values are given in Table [2J In order to assess the 
errors in the parameters we have measured the spectrum of the residuals from the second- 
derivative fit and conducted Monte Carlo simulations by randomly varying the phases of the 
spectral components. Quoted uncertainties correspond to the standard deviations around 
the mean values. We obtain i) = 2.8 ± 0.6 x 10~ 24 Hz s~ 2 , which implies that the present 
braking index of this pulsar is very low, being n = 0.9 ± 0.2. 

Many young pulsars such as the Vela pulsar (PSR B0833-45), PSR B180 0-21 and PSR 
B182 3— 13 have large positive values v associated with recovery from glitches ( iLyne fc Smith 
20061 ). If there has been any recent glitch activity in PSR J1734— 3333, the present value of 
v would very like ly be contaminate d by any such relaxation and the long-term value would 
be even smaller ( jHobbs et al.l 120 101 ). indicating that the actual long-term braking index of 
PSR J1734-3333 must be less than or equal to 0.9 ± 0.2. 



3. Discussion 

This braking index measurement indicates a movement in the P P diagram with a slope 
2 — n of at least 0.9. The inset diagram in Fig. [T] shows how, superposed on timing noise, 
there is a clear and systematic motion over nearly 14 years towards the top-right region 
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of the diagram, where magnetars are found. If this behavior is sustained, we can estimate 
the time it will take for the pulsar to move from its present position in the diagram to the 
magnetar region. Using the power law that defines the braking index (Eq. [JJ), an object 
relocates in Fig. [TJfrom (-Pi, -Pi) to (P2, P2), where P2 = P\{P2 / P\) 2 ~ n , in a time given by 



AT 
AT 



2r e 
n-1 



n-1 



2r c ln ( g 



1). 



(2) 



in 



where r c = P\j1P\ is the present characteristic age. The braking index of a neutron star is 
determined by the mechanism which provides the slow-down torque, due either to electro- 
magn etic radiation or particle flow, and any variation in the moment of inerti a or magnetic 
field ( Manchester et al.l 1 19851 ; iBlandford fc Romanil Il988l ; iLyne fc Smith! 120061 ) . If these un- 
derlying physical processes remain unchanged and PSR J1734— 3333 continues to evolve with 
a constant braking index of n = 0.9, Eq. [2] implies that it would reach a period of 8 s in only 
29 kyr. It would then be situated in the middle of the magnetars. The characteristic age 
would be 6.7 kyr, less than its present value of 8.1 kyr and much smaller than the combined 
durations of its previous life as a radio pulsar and its subsequent path amongst magnetars. 

Assuming the relation between the neutron star magnetic field and slow-down rate given 
in section 11.11 the position in the P-P diagram indicates the effective dipole magnetic field 
strength. Most quoted magnetic field strengths, including those of magnetars, are based 
upon this calculation. The upward movement on the P-P diagram of PSR J1734— 3333 
then corresponds to an increase of the effective surface neutron-star dipole magnetic field. 
In this basic model (magnetic-dipole radiation of an orthogonal rotator in vacuum), n < 3 
can be o btained by variations of eith er the moment of inertia of the star or the magnetic 
moment (IBlandford fc Romanil 119881 ). While the moment of inertia is indeed expected to 
decrease in young pulsars due to re-shaping caused by the secular spin-down, it is unlikely 
that this process can be sustained over long times. In contrast, formation of cracks in the 
crust are expected, that may c ause glitches or abrupt spin irregula rities which would be 
visible in the timing residuals ( IBaym et al.l Il969l ; lAlpar et al.lll996l ). Nothing like this is 
visible in the data for PSR J1734— 3333 during the 13.5 years of observations. Variations of 
the magnetic moment could be caused by alignment or misalignment of the magnetic dipole 
with respect to the spin axis or by increase or decrease of the magnetic field strength. To 
obtain n ~ 0.9 the magnetic axis would have to migrate away from the rotation axis with a 
time-scale of 10 4 yr. No such changes have ever been measured and the available evidence 
for secular magnetic moment migration in radio pulsars points to alignment (rather than 
misal i gnment) on time-sc ales of 10 6-8 yr (ITauris &: Manchesterlll998t IWeltevrede &: Johnston 



2008 



Young et al.ll2010l ). Hence, the only option in the simple magnetic-dipole radiation 
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scenario to produce n < 3 is surface magnetic field growth. 

The above interpretation almost certainly involves significant simplification. Observa- 



Kramer et al. 


2006 


Lvne et al. 


2010) 



Indeed, out-flowing plasmas could remove enough angular momentum from the star to 
become the dominant spin-down mechanism, thereby giving n ~ 1. In this scenario, 
it has been predicted t hat strong relativistic winds c ould produce n = 1 (jMichel 1969 



Manchester et al. 



198 



5; 



Olausen et al. 



Alvarez fe Carraminanal 12004 ; iBucciantini et al.l 120061 ) . We note 



r|2oioh :ound no evidence for extended emission around PSR 



however, that 

J1734— 3333 in X-ray observations, suggesting the absence of a pulsar wind nebula (PWN). 

The ideal magnetic-dipole radiation model has been modified in several occasions to 
produce more realistic models, commonly consisting of a corotating closed-field-line region 
surrounded by an open-field-line region. These models consider different dipole inclina- 
tion angles and incorporate the effects of a plas ma filled magnetosphere (e.g. lMelatoslll997l ; 
Contopoulos fc Spitkovskyl 120061 ; iLi et al.l 1201 ll ). In these scenarios the spin-down can not 
only be regulated by the strength of the dipolar magnetic fiel d, but importantly also by 
the d ipole inclination angle, the spin rate of the op en field lines (Contopoulos fc Spitkovsky 
20061 ). and the conductivity of the magnetosphere (ILi et al.l 1201 ll ) . In these models, a brak- 
ing index of three is still expected to be the norm, while for a braking index n < 3 it has 
to be assumed that either the co- rotating magnetosphere ends well inside th e light cylinde r 
( IContopoulos fc Spitkovskyll2006l ) or the conductivity is increasing with time ( ILi et al.ll201ll ). 
Both may explain the case of PSR Jl 734— 3333, but we remark that the large majority of 
measured braking indices are not three as expected from these models but are smaller than 
three. In the case studied here, the braking index is much smaller than three, so that we also 
consider the alternative explanation that the dipole surface magnetic field may increase with 
time. It is possible, that a combination of effects contributes to the observed spin-down. 

Our ignorance of the actual spin-down mechanism for PSR J1734— 3333 could be partly 
reduced by considering the properties of other high-P objects. The three known radio emit- 
ting magnet ar s exhibit unique radio emission propert i es, not generally o bserved in normal 
radio pulsars ( IKramer et al.l 120071 : ICamilo et al.ll2008l ; iLevin et al.ll2010f ). which have been 
attributed to the presence of strong magnetic fields. Another object, PSR J1846— 0258, is an 
X-ray pulsar that in 2006 showed magnetar-like actiy ity after more than 6 years of normal 
rotation-powered pulsar behavior ( iGavriil et al.l 120081 ) . The effective dipolar magnetic field 
of this pulsar is amongst the highest among normal pulsars and its radiative properties dur- 
ing the 2006 outburst, as with the three radio magnetars, have been related to the effects of a 
very high surface magnetic field. By studying an unusual braking index decrease observed af- 
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ter the 2006 event, iLivingstone et al.l (120111 ) concluded that the effects of parti cle wind losses 



(if pre sent) are not significant in the rotation of PSR J1846— 0258. Finally, lHarding et al. 



( 119991 ) showed (in the case of SGR 1806—20) that even though its rapid slowdown may be 
caused by stellar winds, a high magnetic field may still be necessary to power the typical 
magnetar activity. 

Thus, provided that the singular properties of high-P neutron stars, including mag- 
netars, are effectively driven and determined by very strong magnetic fields, the move- 
ment of PSR J1734— 3333 on the P P diagram could be simply understood as a pas- 
sage to a new magnetar life caused by surface magnetic field growth. The surface mag- 
netic field of a neutron star could increase due to outward diffusion of a stronger internal 
field (IMuslimov fc Pagdll996l ; iGeppert et al.lll999l ; iHol l201ll ). Such a field could have been 
buried by hyper critical accretion of material falling back just after the supernova explosion 
(jChevalierlll989l ) . In the case of ohmic diffusion, this rate will depend mainly on the depth of 
the particular submersion, with calculations suggest ing that the inner field will start emerg 



ing w ith typical timescales between 10 2 to 10 6 yr (jMuslimov fc Pagd 1 19961 ; IGeppert et al. 



19991 ). This is consistent with the time required for PSR J1734— 3333 to reach full magne- 
tar spin properties. It may be that the only modification required to conventional models 
of magnetar formation is a possible delay of perhaps 10-100 kyr in the emergence of the 
magnetic field. It is possible that PSR J1734— 3333 has already undergone some transitory 
magnetar activity and by the time it reaches full magnetar rotational properties it may have 
developed the full radiative properties of a magnetar. 



Conclusions 



The spin evolution of PSR J1734— 3333 during the last 13.5 years is consistent with 
a braking index n = 0.9 ± 0.2, which corresponds to a systematic movement in the P-P 
diagram with a slope of at least +0.9, arriving within the region of the magnetars after 
about 30 kyr. We discussed that secular evolution of the magnetic-dipole inclination and 
systematic variations of the moment of inertia are unlikely to dominate pulsar spin evolu- 
tion. Also, X-ray obs ervations of PSR J1734— 3333 show no evidence for any form of PWN 
( lOlausen et al.l 120101 ) which, if present, would be a signature of wind ac tivity around this 



pulsa r that could affect its rotation leading to a braking index n < 3 (e.g. iBucciantini et al. 
20061 ). Most realistic models involving electromagneti c torques predict n ~ 3, although could 



give n < 3 if the corotating magnetosphere is small (IContopoulos &: Spitkovskvl 120061 ) or if 
the magnetosphere's conductivity is increasing with time ( ILi et al.l 1201 ll ). 



We interpret the movement of PSR J1734— 3333 in the P P diagram, together with 
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the properties of other high-P neutron stars, as evidence for a possible new genesis for at 
least some magnetars. Regardless of whether or not this movement corresponds to actual 
magnetic field growth, our interpretation has several satisfacto ry implications: Firstly , it 
has been suggested that the apparent excessively high birthrate (IKeane fc Kramer! 120081 ) of 
all species of neutron stars over core collapse supernovae events can be explained if there is 
evolution between the species, as we suggest here. Secondly, many magnetars will be older 
than their characteristic ages suggest by 10-100 kyr. Therefore, not all magnetars will be 
young enough to be still surrounded by a visible SNR, consistent with the surprisingly small 



number of secure assoc iations between magnetars and SNRs (Fig. [T]) (IGaensler et al.l 12001 



Allen &: Horvathl 12004 ) . Thirdly, the same physical processes that would be responsible for 



the development of the magnetar properties may be responsible for the effective dipolar 
magnetic field growth seen in the rotation-powered pulsars with measured values of braking 
index, which all have values of n < 3 (Fig.fU Tabled]). Finally, recent observations of Central 
Compact Objects (CCOs) in SNRs suggest that a significant number of ne utron stars are 
born with long periods and low period derivatives ( jHalpern fc Gottheh1l2010l ). Hence, young 
pulsars are not only found at the top-left of the diagram but also towards low P values, 
populating the left hand side of the main bulk of pulsars and probably being progenitors 
for most normal pulsars. This may imply that the high-P young pulsar population are not 
the progenitors of most pulsars, making feasible their possible evolution into magnetars. In 
summary, we suggest that, with a range of internal magnetic fields, submersion depths and 
initial spin periods, we could explain the e xistence of m agnetars, CCOs and the whole radio 
pulsar population as one single family (cf. lKaspill2010l ). 
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Fig. 1. — P-P diagram for all known magnetars and young pulsars having P > 1.65 x 10~ 15 . 
Lines of constant magnetic field are dashed and lines of constant characteristic age are 
dotted, with a slope of 1. Open arrows indicate the motion of those pulsars with published 
values of braking index (Tabled]). Each arrow represents the projected motion of the pulsar 
during the next 10,000yr (excepting the one for PSR J0537— 6910, for which only 2,000yr 
was used), assuming that it evolves with a constant braking index n. Equal logarithmic scales 
are chosen, so that a pulsar will move with a slope of 2-n. The motion of PSR J1734— 3333 
is represented by a closed arrow. Objects which have robust proposed associations with 
supernova remnants or pulsar wind nebulae are identified with a surrounding large square, 
with smaller squares for less convincing associations. White triangles are used for the radio 
emitti ng; magnetars. The P value of the SGR with the lowest P corresponds to an upper 
limit (IRea et al.ll2010i ). Inset: the pulsar's motion in P P space over a 13.5-year period using 
equal logarithmic scales, but magnified by a factor of 4,000. Major ticks are separated by 
5 x 10 -4 s on the horizontal axis, and by 3 x 10~ 4 on the vertical axis. Each data point in the 
inset is the result of a fit of P and P to a 1500-day interval set of times of arrival. The error 
bars are the standard deviations (See section [2]). Most of the information was taken from the 
ATNF pulsar catalogue (version 1.39, |http://www.atnf.csiro.au/research/pulsar/psrcat/ ), 
or the McGill SGR/AXP Online Catalog as it was in November 2010 
( http: / / www.physics.mcgill.ca/ ~pulsar/magnetar/main.html[ ) . 
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Fig. 2. — Timing residuals of PSR J1734— 3333 over a 13. 5- year interval. The timing residuals 
are the difference between the observed times of arrival of pulses and those predicted by a 
simple slow-down model of frequency and first derivative v. The curve shows the slow-down 
model after including the fitted value of frequency second derivative of 2.8 ± 0.6 x 10~ 24 Hz 
s -2 (Table [2]), corresponding to a braking index of n = 0.9 ± 0.2. 
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Table 1. Published braking indices. 



Pulsar 


n 


Reference 


B0531+21 


2.51(1) 


Lvne et al. ( 19931 


J0537-6910 


-1.5 


Middleditch et al. (2006) 


B0540-69 


2.140(9) 


Livingstone et al. (2007) 


B0833-45 


1.4(2) 


Lvne et al. (1996) 


J1119-6127 


2.91(5) 


Weltevrede et al. (2011) 


B1509-58 


2.839(1) 


Livingstone et al. (2007) 


J1846-0258 


2.65(1) 


Livingstone et al. (2007) 


J1734-3333 


0.9(2) 


this work 



Note. — The uncertainty in the last shown digit is quoted in parenthesis. While all other 
values come from phase- coherent timing, the values for PSR B0833— 45 and PSR J0537— 6910 
were obtained by studying the evolution of the spin-down rate after several glitches. The 
braking index for this last pulsar was calculated from a rough estimate of the systematic 
variation during more than 8 yr of the frequency derivative data. We note that the braking 
index of PSR 1 846—0258 was found to b e 2.2 ±0.1 when using data following the 2006 event 
(see section 2) (jLivingstone et al.ll201ll ). 
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Table 2. Observed and derived parameters for PSR J1734— 3333 obtained from fits to the 

pulse times-of-arrival. 



Parameter 



Value 



-15 



Hz s 
Hz s 



- 1 ) 

- 2 ) 



RAJ 
DECJ 

v (Hz) 
v (10" 
v (10~ 24 
P(s) 

P (io- 15 ) 

P (10- 24 s 
Timing Epoch (MJD) 
Data span (MJD) 
DM (cm~ 3 pc) 
Suoo (mJy) 
W 50 (ms) 



17:34:26.9(2) 
-33:33:20(10) 
0.855182765(3) 
-1667.02(3) 
2.8(6) 
1.169340684(4) 
2279.41(4) 
5.0(8) 
53145 
50686-55602 
578(9) 
0.5 
500 



Distance (kpc) 7.4 

Characteristic Age (kyr) 8.1 

Surface magnetic field (TG) 52 

Braking Index, n 0.9(2) 



Note. - Standard errors are given in parenthesis in units of the last quoted digit. See 
section [2] for more details. 



